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Abstract 



In the presence of noncompact symmetry, the stability of relative equilibria 
under momentum-preserving perturbations does not generally imply robust sta- 
bility under momentum-changing perturbations. For axisymmetric relative equi- 
libria of Hamiltonian systems with Euclidean symmetry, we investigate different 
^ mechanisms of stability: stability by energy-momentum confinement, KAM, and 

^SJ Nekhoroshev stability, and we explain the transitions between these. We apply 

our results to the Kirchhoff model for the motion of an axisymmetric underwater 
vehicle, and we numerically study dissipation induced instability of KAM stable 
relative equilibria for this system. 

1 Introduction 

O 

> 



Relative equilibria of Hamiltonian systems with symmetry are special solutions which 
are equilibria of the symmetry reduced dynamics. Relative equilibria are called stable 
if they are stable equilibria for the symmetry reduced dynamics. In other words, in 
a Hamiltonian system with symmetry group G, a relative equilibrium is called stable 
(more precisely ^-stable [13,14,16]) if every time-orbit close to the relative equilibrium 
stays close the the fj-orbit of the relative equilibrium. Since the fj-orbit of the relative 
equilibrium usually strictly contains its time-orbit, tj-stability is usually weaker than 
orbital stability. 

The stability of relative equilibria is delicate for Hamiltonian systems with non- 
compact symmetry. We have established |14j exactly why it is erroneous in general to 
conclude stability of relative equilibria under momentum-changing perturbations from 
stability under momentum-preserving perturbations. In the presence of noncompact 
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symmetry, there is a gap between these. In that gap, energy- momentum confine- 
ment fails, meaning that stability under momentum-changing perturbations cannot 
be established by energy-momentum Lyapunov functions. 

Here we investigate the stability of axisymmetric relative equilibria of Hamiltonian 
systems which admit the Euclidean symmetry Q = S0(2) k M'^ x S0(2). Our work 
bears on some results of Leonard and Marsden [7] concerning a class of relative equilib- 
ria in the 12-dimensional Kirchhoff model for an axially symmetric underwater vehicle, 
in which the vehicle falls and spins. They derive a condition for energy-momentum 
confinement under arbitrary perturbations ([7], Theorem 4.4), and they calculate that, 
at the boundary of this stability region, a Hamiltonian Hopf bifurcation occurs. We 
show that energy-momentum confinement under arbitrary perturbations actually oc- 
curs in a smaller region, but that, in the intervening gap, stability can be established 
by KAM methods. We also show that the coincidence of the gap boundary and the 
Hamiltonian-Hopf bifurcation is to be generally expected when the phase space is 
1 2-dimensional . 

In the underwater vehicle system, we provide numerical evidence that the KAM 
stability in the gap is destroyed by small dissipation, whereas the stability by energy- 
momentum confinement is preserved. In this system, the transition from energy- 
momentum region to the gap is spin independent, whereas the transition from the gap 
to spectral instability, at which the Hopf bifurcation occurs, does depend on spin. All 
the axisymmetric relative equilibria that are spin-stabilized are in the gap, and their 
KAM stability is destroyed by small momentum-preserving dissipation. The implica- 
tion for gyroscopically stabilized devices is startling. In the presence of noncompact 
symmetry, robust stability may not be achievable by the use of spin: dissipation induced 
loss of stability of the relative equilibrium will occur even in absence of dissipation of 
spin. When the symmetry is noncompact, a general understanding of the kinds of 
stability which generically occur, and the transitions between them, is necessary for 
the determination of robust stability criteria. 

This paper is structured as follows: We begin in Section [2] with a description of ax- 
isymmetric relative equilibria of Hamiltonian systems with Euclidean symmetry. We 
introduce coordinates related to the reduction of the system by its largest abelian sub- 
group /C = M'^ X S0(2) of G- We obtain a family of Hamiltonian systems parametrized 
by the corresponding momenta. This lays plain the essential difficulty, because it shows 
that perturbations to arbitrary momentum are S0(2) symmetry-breaking, for an addi- 
tional SO (2) symmetry present in the system which does not commute with the abelian 
subgroup /C. The symmetry-breaking occurs because the additional S0(2) symmetry 
acts on the whole family of Hamiltonian systems, including the parameterizing mo- 
menta, and acts on one single Hamiltonian system of the parameterized family only 
for certain vertical momenta; only those Hamiltonian systems are S0(2)-symmetric. 
In Section [3] we study the stability of axisymmetric relative equilibria. In general, 
we establish the presence of the gap and discuss different mechanisms of stability, 
and for Hamiltonian systems with 12-dimensional phase space, we establish the Hopf 
eigenvalue collision at the gap boundary. In Section [4] we begin with a brief summary 
of the Kirchhoff model. We recover the stability criteria of J^ for the falling, spin- 
ning relative equilibria, and prove stability of the relative equilibria within the gap by 
verifying the Moser twist condition for the corresponding equilibrium on the reduced 
space. In Section [5] we numerically demonstrate that addition of small dissipation 
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distinguishes energy-momentum from KAM stability, by observing that stabihty is 
maintained in the former and destroyed in the latter by the addition of momentum 
conserving dissipation. 
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2 Dynamics near axisymmetric relative equilibria 

In this section we introduce axisymmetric relative equilibria of Hamiltonian systems, 
and introduce suitable coordinates near them for use in the stability analysis in Sec- 
tion [3] We consider a general context that includes the underwater vehicle example 
in Section m Let 

X = fnix) (2.1) 

be a Hamiltonian system defined by an energy : ^ M on a connected symplectic 
manifold (Al,a;), i.e. 

uj{x)[fH{x),w) = DH{x)w for all x e M,w e %M. 

Let ei, 62, ea denote the unit vectors along the x, y, z axis of M.^, and for any w e M'^ 
denote 



so that, for example. 
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Define the group 

g = SE(2) X M X S0(2) = S0(2) x x S0(2), 

where x stands for a semidirect product in that the multiplication for g = ((f), a, 9) G 
S0(2) X M3 X S0(2) is 

gi92 = {(t)i,ai,9i){(l)2,a2,02) = Ui + 02, «! -I- exp(e3(/)i)a2, 6*1 + 92). 
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Assume G acts symplectically on M, i.e. u> is ^-invariant, and suppose that H is G- 
invariant. This imphes that the vector field Jh commutes with G, as does its flow, i.e. 
Jh and its flow are C/-equi variant. Also, assume that G acts properly and faithfully. 

In the following, S0(2)/j is the copy of S0(2) in G which lies in the center of G, 
so that e e S0(2)i^, and S0(2)i is such that (p e 80(2)^. Here "i?" stands for 
right and "L" for left. In the application to axisymmetric underwater vehicles, see 
Section [ijbelow, the S0(2)l K is related to spatial isotropy, it corresponds to a left 
multiplication, and it translates and rotates the body in space. The S0(2)fj action is 
related to a material symmetry, it corresponds to a right multiplication, and it spins 
the body around a symmetry axis. 

By Noether's theorem there are dimC/ — 5 locally defined conserved quantities 



of (2.11. We assume these exist globally and organize them in the momentum map 
J: M ^ Q* such that uj{x){^x, w) = DJ^{x)w for all x G M, w G %M, £, € g. Here 
— so{2)l k X so(2)fl; = TidG is the Lie algebra of G- Denote the components of 
J by 

J(x) = (J*, J^ J'')(x), J*(x), J^(x) eso(2)*, J"(a;)e(M3)* (2.3) 

and the elements of g* by 

M = (M^ /) eg* - 50(2)1 X K3 X so(2):^. 

In the underwater vehicle example of Section |4] the component J" is the linear mo- 
mentum, J"* the angular momentum, and the momentum of the spin of the body. 
As in that system, we will assume that the momentum mapping J transforms by the 
coadjoint action on g* , i.e. [8 

J{gx) = (Ad*)"^J(x) for any g eG and x e M. 

Here Ad* is defined by the condition (Ad* = n{Adg 0, fJ- & g* , £, & g, g £ G , and 
Adg^ — g£,g~^- A standard computation (see, for example, [5]) gives 

J*((0,a,0):E) = J^{x) ~ (exp(e30)a)2J"^(a;) + (exp(e30)a)i J^^a;), 
J"' {{(p, a, e)x) = J"' {x) cos (/) - J'^^ {x) sin 

J""'- ((0, a, 9)x) = J"^ {x) cos (t> + J"^ (x) sin (2.4) 

r-'{{<ly,a,d)x)^r-{x), 

J\{(f>,a,0)x) ^ J^x). 

2.1 Axisymmetric relative equilibria of Euclidean group ac- 
tions 



Let Xe be a relative equilibrium of (2.11, i.e. let 

= (Cf , C, e S0(2) ® E3 e 50(2) = 0, 

and suppose exp(^ei)a^e is a solution of the differential equations x — fuix). We 
assume Xg is axisymmetric, meaning that its isotropy subgroup Gx^ — {g ^ Q • 9^e = 
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Xe} is the diagonal subgroup 



SO(2)i5 = {(0,0,0) e S0(2) K X S0(2)} C g. 

As is general for Hamiltonian systems with symmetry, if rje is in the isotropy algebra 
fla,^ then (^e + rie)xe = ^eXe = fnixe) because r]eXe = 0, hence exp(i(^e + rie))xe is 
also a solution of the Hamiltonian system. Thus is only determined up to addition 
of elements in 

0,^ = so(2)d = { (e^ 0, e'^) e so(2) K M3 X so(2) } C g, 

and so we can choose ^1" = 0. In the underwater vehicle, Xe is an equilibrium in a 
frame that co- moves with the action of exp(^e^), 1^"^°* = ~ = is the angular 
velocity, and is the translational velocity. The momentum value = J{xe) of 
the relative equilibrium is of the form /ig = (/if , /x" , /^g) with /i" || 63, because it is 
fixed under the action of S0{2)o that occurs in the transformation rule (2.4). The 
momentum map jSO(2)d q£ ^j-^g action of the symmetry group S0(2)£) is 

It follows from = that X Q IS cl relative equilibrium for the abelian symmetry 
group 

JC = M!^ y. S0(2)k = { (0, a, 0) G S0(2) k x S0(2) } C (2.5) 

and Xe becomes an equilibrium after reduction by this group. The stability of this 
equilibrium implies stability of the relative equilibrium, see Section |2.2| below. We 
will study the stability of such an axisymmetric relative equilibrium, which in the 
underwater vehicle example of Section [4] corresponds to a vehicle spinning about its 
symmetry axis with angular velocity 0;'^°*, and translating along its symmetry axis 
with translational velocity . 

2.2 Reduction by /C = x 80(2)^ 

The following theorem provides the coordinates which we require for the subsequent 
stability analysis. 

Theorem 2.1. In a Q-invariant neighborhood oj an axisymmetric relative equilibrium 
Xe of (2.1 1 there are coordinates x = {a,6,h'"',i'^,w) such that 



(a) J" = J = V , and the differential equations (2.1 1 are 

ii = T)^.H{v,w), e = T)^bH{v,w), v'^^Q, v^^Q, w = IT)^H{v,w), (2.6) 

where v ~ (v°'^v^), J is the standard symplectic structure m,atrix on the linear 
symplectic space W — KdiinA4-8^ ^^^^ ^ ^ 
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(b) The Hamiltonian H^VjW) is invariant under the action o/SO(2)/), which takes 
the form 



v°- exp((^e3)z^°, , w — > Rc^iy, w), 

i.e. the action on w generally depends on v and this equation defines R^. 
(c) The coordinates of the axisymmetric relative equilibrium x^. are 

a=0, 6' = 0, V^Ve^ {j''{Xe),j\xe)), W = Q, 



and w — is an equilibrium of the w-equation of (2.6) at v = v^. 

The a — Q, = Q plane is locally a slice at Xe through the action of /C = M"^ x 80(2)/? 
and models the Poisson reduced space Ai/JC. The symplectic leaves of this Poisson 
space are given by fixing v. At fixed v, the system 

w^mn,H{v,w) (2.7) 

is the (Marsden-Weinstein) symplectic reduced system [5] corresponding to the abelian 
subgroup /C at momentum u, obtained simply by "ignoring cyclic coordinates" . The 
additional 1-dimensional symmetry is expected because /C is a codimension 1 subgroup 
of Q. Most of the proof follows from the general theory developed in [15], but since 
the symmetry group K, is abelian there is also the following elementary proof. 



Proof of Theorem 2.1 The four momenta J^, J 2, Ja, J Poisson commute since K, is 
abelian, and their derivatives are linearly independent at Xe since the isotropy group 
of Xe is not bigger than SO(2)i3. Define v"' — J°,i^^ ~ J^, and using the Darboux 
theorem |18j . choose four conjugate functions, i.e. four functions ai, 02, 03, 6, such that 

Find functions qi, . . . ,qk,Pi, . ■ ■ ,Pk, k = ^(dimAI — 8), that Poisson commute with 
the original 8 and satisfy 

= {%,9j} = 0, {p^,PJ} = 0. 

Set w — {q,p), and translate so that Xe is at w = 0. The group /C acts by addition on 
the conjugate coordinates 01,02,03,^ since this action is generated by the momenta 
i/O-i ^ 1/^2 ^ jyO-a ^ j^e ^ So H does not depend on the coordinates ai,a2,a^,0, and we have 
the v — {v°'^v^) parametrized canonical system on the linear symplectic space W 
given by the system (2.7 1 of (2.6 1. The other equations of (2.6 1 follow immediately. 

The symmetry group S0(2)£) acts on the variables {Wjv) independently of a, 9. 
Indeed, if / is a function which does not depend on a, 6, if x is one of ai, 02, 03, and 
the conjugate of a; is a;, then 

^{j'' + J^/} = {a;,{j'^ + J^/}} 

= -{ J'^ + J^ {/, x}} - {/, {i, J' + J^}. (2.8) 
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The first term of (2.8 1 is zero because 



5/ 

dx 



For the second term of (2.8 1, note that x is one of v"'^ , v""^ , v""^ , v , which are mo- 
menta of /C. Since /C is a normal subgroup of Q, its momenta are a Poisson ideal of 



the momenta J^, C € fl. This can also be seen directly from (2.4 1 because the Poisson 
brackets of J*^ with the other components of the momentum map J are found by 
differentiation in at = of the right hand sides of (2.4). The results are linear 

- J"^} is a linear combina- 



combinations of J and J° , i 



1,2,3. In particular, {x^J 

is zero for the same reason as the 



, so the second term of (2 



tion of v°-^ , v°-' 

first. Thus there is an action of S0(2)£i on the variables {v, w). Since J is equivariant 
and J — {J'^,v), the resulting action on the v variables is independent of w and 
equal to the coadjoint action of (6', 0,6*) on v, which by (2.4 1 is a rotation of {^1,1^2) 
by 6. The coordinates above can be restricted to an S0(2)£i invariant neighborhood 
because Xe is fixed by the action of this compact group and are therefore coordinates 
in a C/-invariant neighbourhood oi Qx^. □ 



Definition 2.2. A momentum v G 
coadjoint action o/SO(2)£), i.e. if v"^ 



63 e 



so(2))' 

10)3 



IS 



vertical if it is fixed by the 



For vertical v, the system (2.7 1 inherits an S0(2)£) symmetry from the full phase 
space M. It does not have this additional S0{2)d symmetry for nonvertical v. Hence 
a perturbation of the S0(2)£)-symmetric Hamiltonian system (2. 7 1 at a vertical mo- 



mentum value v to a, non- vertical momentum value is an S0(2)£)-symmetry breaking 
perturbation. 



Proposition 2.3. For vertical momenta v, 
with momentum map 



the S0{2)]j action on W is symplectic 



'w 



0=0, 9=0 



In particular, at vertical momentum u, the system (2.7 1 is S0(2)£i symmetric and 
conserves Jw- 



2.1 



IS 



Proof. The action of S0(2)o in the coordinates {a,6,v°',i' ,w) of Theorem 
symplectic since it is the action of a subgroup of Q, which acts symplectically by 
assumption. It follows that S0{2)d acts symplectically on W, since the symplectic 
form on W at fixed v is the restriction of the symplectic form on A4. Using the 



symmetry properties (2.4 1 of the momentum map 



= (J* + J')((-a,-0)(a,0,I/^^/^w;)) 
= (J-^ + J^){a, 9, i/", iy^,w) + a^vl - a^v^ 



(2.9) 



so Jw — jSO(2)£5 j£ _ _ g xhus for fixcd vertical i/, Jw is equal to the 
momentum generating the S0{2)u action on W. □ 
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Remark 2.4. Equation (2.9 1 and the differential equation for a in the system (2.6 1 
imphes that the conservation of Jw is not typical for nonvertical momenta. Indeed, if 
Jw is conserved at some fixed nonvertical ly = uq, then adding any function eH{v, w) 
to H results in 

^(a2^^o^ - "li'os) = e I ^^01 " ^^^02 1 • (2.10) 




We can choose an S0(2)c invariant function H such that the right side of (2.101 is 
not everywhere zero for e arbitrarily small. 

Remark 2.5. Locally, the symmetry reduced space takes the form {1^,10) t* W 
where t is the Lie algebra of /C. The space W is called the symplectic normal space 
at Xe with respect to the symmetry group /C, see |15j . since it is transversal to the 
group orbit )Cxe at Xe and is a symplectic space. Moreover it is the largest symplectic 
subspace of the normal space J\f to the group orbit Qx^ at Xe- In the notation of [15], 
we have A/q = i* , Ni = W. An application of the results of [T5] would mean a 
reduction by the group H = SO (2) ad x where 

S0{2) AD = {(0,0,-0) £ S0(2) kR^x S0(2)} 

is the antidiagonal embedding of S0(2) in S0(2) l x 80(2)^. But the resulting Poisson 
structure on A/o = f)* — se(2)* ©K* (where f) is the Lie algebra of H) is nontrivial since 
the symmetry group H is nonabelian, so that M ~ A/q ® A/i. In particular the sym- 
plectic leaves have nonconstant dimension at the axisymmetric relative equilibrium. 
For this reduction the conserved quantity jSO(2)i3 ^ Poisson momentum map for the 
S0(2) D action on Ao ® A/i, and this momentum map is a conserved quantity for the 
reduced system on Ao ® Afi ■ But the nontrivial Poisson structure on Ao would make 
the stability analysis more difficult. We would have to employ the general methods 
developed in |14j to study stability by energy-momentum confinement and we would 
have to desingularize the Poisson structure using blow up methods as in |13j to study 
the KAM stability of the axisymmetric relative equilibrium Xe- 



Remark 2.6. In Theorem 2.1^ we could also construct an S0(2)o invariant Witt 
decomposition of the tangent space T^^M and then use the equivariant Darboux 
theorem to obtain a model IC x t* x W oi a, tj-invariant neighborhood of Qxe- In this 
construction the S0(2)£) action on W would be linear and ;/-independent, see [15J. 
Note however that in this case i* is not the annihilator of g^;^ = 50(2)^1, so the situation 
here is not exactly the setting of |15j . 



3 Stability of axisymmetric relative equilibria 



In this section we use the coordinates from Theorem |2.1| to study the stability of the 
axisymmetric relative equilibrium x^- We show that dcfinitcness of the Hessian of 



H{i'e,w) at the equilibrium w = of (2.7\ implies stability of Xe under arbitrary 



perturbations. When this Hessian is not definite, for momentum-preserving pertur- 
bations, a constant times the momentum Jyy can be added to the energy to obtain 



a Lyapunov function of the system (2.7). However, for nonvertical linear momentum. 
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Jw is not a conserved quantity of (2.7 1, so it cannot be used in a Lyapunov func- 



tion to establish stability for perturbations to nonvertical linear momentum. The gap 
(Sections [l] and 3.1 ) is defined to be those relative equilibria for which a constant 
times the momentum J^y must be added to obtain a Lyapunov function. If Xe is in 
the gap then the equilibrium w = of (2.7 1 has imaginary spectrum and indefinite 



Hessian, and is amenable to classical KAM and Nekhoroshev results (Section 3.11 to 
obtain stability. 



In the case of a 12-dimensional phase space A4, the system (2.7 1 has 2 degrees 



of freedom, and Propositions 3.9 3.11 below show that formal stability (Section 3.11 
of the axisymmetric relative equilibrium fails exactly at an eigenvalue collision. 
Moreover we show that this is exactly where a Hamiltonian Hopf bifurcation typically 
occurs and a loss of linear stability is expected. Thus, in the case of a 12-dimensional 
phase space, formal stability and linearized stability of axisymmetric relative equilibria 
are typically equivalent, explaining the observation of this by Leonard and Marsden 
in the Kirchhoff model of axisymmetric underwater vehicles [7] . 



3.1 Nonlinear Stability 

The relative equilibrium Xg is called A stable for some subset A oi Q, if initial data 
starting sufficiently close to Xg stay arbitrarily close to Axe for all times. It is formally 
stable if TP'{H + \J\Y){xe)\w is definite for some A G M. Recall from (2.5) that 
/C = X S0(2)fl. 



Proposition 3.1. The axisymmetric relative equilibrium Xe of (2.1 1 is JC-stable if the 
Hessian T)^H{iyf,,0) is definite. 



Proof. If the Hessian D^i/(:/e,0) of the Hamiltonian system (2.7 1 is definite then, 
using energy as a Lyapunov function, the equilibrium w = of {2.7 \ is stable under 
perturbations in w, i.e. within = i^g. The equilibrium w — persists to nearby 
momentum values and the corresponding equilibria of (|2.7| at those momentum 



values are also stable. Since the momentum i/{t) is a conserved quantity of (2.6) we 
conclude that, for initial data close to i/ = v^tW = 0, the solution (i/(t), w(i)) of (2.6) 
stays close to (i^ejO) for all times. Since the only remaining variables are a, 0, and 
these correspond to the action of /C, the proof is complete. □ 

Remark 3.2. Better stability criteria cannot be obtained by relaxing to the weaker 
stabihty modulo the whole group G, instead of just the subgroup /C, because the 
isotropy of Xe implies that these two stabilities are equivalent. 

Let j'^ = (J°, J^) denote the momentum map for the symplectic action of K. on 

M. 



Proposition 3.3. The axisymmetric relative equilibrium x^ of (2.1) is K.-stable under 



perturbations that preserve the momentum J of K if it is formally stable. 

Proof. A perturbation which preserves j'^ does not change v and therefore is a per- 
turbation of the stable equilibrium w = Q of the system (2.7). Since Jw is a conserved 



quantity of (2.7) at v — Vg by Proposition 2.3 the equilibrium w = of (2.7) is stable 
if there is some A e M such that T)'^H{Q) + AD^ Jvt'(O) is definite. □ 
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So we see that the presence of the SO (2) d symmetry and its momentum Jw{w) at 
vertical momenta || 63 causes a gap between energy-momentum confinement under 
momentum-preserving and general perturbations of Xe- 

Definition 3.4. Let Ae be the set of axially symmetric relative equilibria. The EM- 
region is the subset C Ae such that D^iJ(z^e, 0) is definite. The gap is the subset 
_4gap ^ ffifj^f formally stable but not in the EM-region. 

Remark 3.5. If the symmetry group Q is compact then the gap is absent since formal 
stability implies tj-stability, see e.g. [T^IIIJ- 

We will say that the axisymmetric relative equilibrium Xe is nondegenerate if the 
Hessian D^_ff(z^e,0) of the equilibrium it; = of (2.7l is invertible. Suppose that 
Xe G Af^^ is nondegenerate. Then w = is a stable equilibrium of (2.7 1, so the 
linearization JD^_ff(i^e,0) of (2.7 1 has nonzero purely imaginary eigenvalues, i.e. it 
is elliptic. We call the relative equilibrium Xe of (2.1) elliptic if the corresponding 
equilibrium it; = of ( |2.7| is elliptic. Depending on the dimension of the phase 
space M, difl^erent scenarios are possible. 



(a) If dim = 10 then ( 2.7 ) is a 1 degree of freedom Hamiltonian system. Equilibria 



of such Hamiltonian systems have definite Hessians if and only if they are elliptic. 
So, the Hessian D^7J(i/e,0) is definite, contradicting Xe G Af^^. In particular, 
at this dimension there are no nondegenerate relative equilibria in the gap. 



(b) If dimAl ~ 12 then, for each v, (2.7 1 is a 2 degree of freedom Hamiltonian 



system. The equilibrium w — persists as elliptic equilibrium of (2.7 1 to nearby 



v. Elliptic equilibria of such Hamiltonian systems are expected to be KAM 
stable [2], and so the axisymmetric relative equilibrium Xe is expected to be 
/C-stable. 

(c) If dim Al > 14 then the equilibrium w — persists as in case (b), but for each 
ly the system (2.7 1 has at least three degrees of freedom. The equilibrium w = 



of (2.7 1 is therefore not expected to be stable because it will exhibit Arnold 
diffusion. However, Nekhoroshev stability |5j [IT] is expected, which implies 
/C stability of the axisymmetric relative equilibrium Xe over exponentially long 
times. 

Remark 3.6. In case (b), in order to prove KAM-stability of the axisymmetric relative 
equilibrium Xe, one can verify the twist condition (see [9 ) for the equilibrium w — of 
the 4-dimensional S0(2) symmetric Hamiltonian system (2.7 1 at = Ve. If the twist 
condition holds ai v — i>e, then it holds for the equilibria persisting to nearby v and Xe 
is therefore /C-stable. The additional S0(2) symmetry of (2.7) at v = Ve can be used 



to express the reduced Hamiltonian of (2.7) in terms of S0(2)/)-invariant functions 



This simplifies the verification of the twist condition, see Section [4~4| But, as already 
emphasized, the conserved quantity Jw of (2.7) at v = Ve cannot be used because it 



is only conserved for vertical v. The system (2.7) for which the twist condition needs 



to be verified, can be obtained in a concrete example by choosing any realization of 
the M.^ X S0(2)j^ Marsden-Weinstein reduced space. 
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Remark 3.7. If Xe € Af^'^ then, as already noted above, Xe is a spectrally stable equi- 



librium with indefinite Hessian of the Hamiltonian system (2.7 1. As shown in O H], 
adding arbitrarily small dissipation to such Hamiltonian systems results in spectral 
instability, because the negative eigendirection of the Hessian forces instability in the 
presence of decreasing energy. Consequently, adding a small S0(2)d invariant dissipa- 
tion for each vertical v, and extending that to nonvertical v, will result in dissipation 
induced instability at each persisting equilibrium of the Hamiltonian systems param- 
eterized by v. If the dissipation preserves also the S0(2)£) momentum Jw, then it 
preserves all momenta resulting from the symmetries of the original system on A^. In 
this case, and if D'^{H + XJw){xe) is positive definite for some A, then the persisting 



equilibria of (2.7) remain spectrally stable for nearby vertical v while, for nearby non- 
vertical ly, there is again spectral instability. In Section [s] we verify this by numerical 
simulation of an x 80(2)^ reduction of the underwater vehicle system. 

3.2 Hamiltonian Hopf bifurcation 

In this section we assume that the phase space A4 is 12-dimensional, so that the 



system (2.7 1 is 4-dimensional. The linearization of the action of S0{2)]j on W at 
the equilibrium w = is a linear symplectic representation, about which we need a 
few elementary facts [10 . Consider a symplectic representation of S0(2) on a 2d- 
dimensional linear symplectic space W. These representations are classified by tuples 
of integers n-i < n2 < ■ ■ ■ < Ud- We say that the action is of type (ni, 712, ... , rid) if 
there is a linear splitting of W into the sum of two 2-dimensional invariant subspaces 
on which the action is isomorphic to z i— > e*"^^z on the vector space = C with 
its standard symplectic structure. With respect to such a splitting, the associated 
momentum mapping is 

^"-ikip H 1- ^'T'dkdP, 

and particularly, the quadratic momentum map for type (ni,n2, . . . ,nd) actions has 
both positive and negative definite eigendirections if both positive and negative rij 
occur in its type. 

Suppose W is a linear symplectic space and A: W — > W is infinitcsimally symplec- 
tic. Let ±iuj be a simple complex conjugate pair of purely imaginary eigenvalues of A. 
The corresponding (2-dimensional) real eigenspace is invariant for the symplectic fiow 
exp(At), which defines a symplectic representation of S0(2) of either type -1-1 or type 
— 1. The Krein sign of ztiui is the sign of its type. We adopt the convention that the 
frequency w has the same sign as the Krein sign of the eigenvalue ±ia;, so that the 
corresponding quadratic Hamiltonian on the real eigenspace, with respect to a basis 
with the standard symplectic structure, is +^uj{q'^ + P^)- Since the frequencies are 
signed, resonances between two frequencies are also signed e.g. there can be —1 : 1 
and 1 : 1 resonances, and these are distinct cases. 

Remark 3.8. Actions of S0(2) of type ( — 1,1) occur frequently in dimension 4, and 



this case occurs in the underwater vehicle example, see (4.251. This case can be 



realized as q t—^ R,pq, p 1— > R^iP where i?^ is the action of S0(2) by counterclockwise 
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rotations (see the proof of Proposition 3.11 1. This is the standard S0(2) action on 



p4 



obtained by lifting the standard action of SO (2) on = {q} to the cotangent bundle 

Proposition 3.9. Let be an axisymmetric relative equilibrium of the S0(2) K 
X S0(2) symmetric Hamiltonian system, (2.1 ) on a 12-dimensional phase space A4. 



Assume that x^ is semisimple and elliptic, and let the action o/SO(2) be of type rij on 
the real±iujj eigenspace, where uji,uj2 are the two normal frequencies of the equilibrium 
w — Q. Then Xe is formally stable if and only if there is not an ni : n2-resonance 
between lui,uj2, ie. if and only if ujin2 — uj2ni ^ 0. 



Proposition |3.9| implies that, in 12 dimensions, elliptic axisymmetric relative equi- 
libria are typically formally stable because resonance is atypical in the linearization 
of the corresponding equilibrium in the reduction by /C. The proof of this Proposi- 
tion follows by applying the following result to the linearized SO (2)13 action at x^, 
restricted to W . 

Proposition 3.10. Let W be a ^-dimensional linear symplectic space and suppose 
S0(2) acts linearly and symplectically on W with momentum map J . Suppose that 
H : W —^ M. is SO (2) invariant and that € W is an elliptic semisimple equilibrium 
where the linearization JD^_ff(0) has eigenvalues iujj, j = 1,2. Suppose that o^i < 
and UJ2 > 0, and let the action of S0(2) be of type nj on the real ±iujj eigenspace. 
Then there is a X such that D^iJ(O) + AD^J(O) is definite (i.e. ^ W is formally 
stable) if and only if 

ujin2 - L02ni ^ 0. (3.1) 

Proof. Split W = into the real eigenspaces of JD^i7(0) to the eigenvalues 
iojj, choosing a basis giving the standard symplectic structure matrix on each factor. 
Then for w = (^1,^2) G ® 

{w, {B^H{0) + XD^J{0))w) = ^(c^i|u.i|2 + ^2^2!') + ^("ikiP + n2\w2\^) 

^ ^ (3.2) 

= -{{uji + A7ii)|wip + (cj2 + An2)|w2n. 

Since the action of S0(2) is nontrivial, one of rti,rt2 is not zero, and 

A 1-^ (^1,^2) + A(ni, 712) 



is a line. To show that (3.2 1 is definite, its is sufficient that this line meets the first or 



third quadrant. But such a line is contained in the second and fourth quadrant only 



if it contains the origin, which is excluded by (3.1 ) 



Conversely, suppose that LUin2 — W2't-i = 0. Then there is a k such that (rii, ^2) 



^(0^1,^2) and (3.2 1 becomes 



(w, {B^H{0) + XD^J{0))w) - i(K + A)(c^i|u.i|2-f c^2h2p), 
and this is not definite for any A because wi < < 072. CH 
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We now consider the case when condition (3.1 1 is violated. 



Proposition 3.11. Assume the setting of Proposition 3.10 again, hut let H{y,w) 
depend on a parameter S K and assume, as before, that the SO (2) action is of type 
(711,712). Assume that 

(i) the equilibrium is elliptic for u <Q, v ^ 0; 



(a) condition (3.11 holds for v <0, ly 0, but is violated at v = Q. 
Then: 

(a) if ni 7^ — 7i2 then the equilibrium is linearly stable for v > Q, v kQ; and 



(b) if Til = —712 BL'^'d if the transversality condition (3.3 I holds then a Hamiltonian 
Hopf bifurcation occurs and the equilibrium becomes spectrally unstable for 
u>0. 

Proof. For part (a), if ni ^ ±712 then the condition LUin2 — 102^1 = Oatw = 0, i^ = 
imphes that the eigenvalues iojj do not collide at = 0. Therefore the symplectic 
eigenvalue theorem pj implies that the eigenvalues are imaginary in a neighborhood 
of v = 0. If 77i = 71,2 then = wi772 — uj2ni = 771(1^1 — L02). Since the S0(2) action 
is assumed to be nontrivial, we have ni ^ 0, and so uji ^ uj2- This contradicts the 
assumption lji < < ^2- 

For part (b), so assuming 77 = 712 = — 77i, there is aequivariant Lagrangian splitting 
W = 1^2 ^^^^ ^Yiai the SO (2) action is z e^^'z on each factor ^ C. 
The standard symplectic structure matrix blocks with respect to this splitting, is the 
identity on the (l,2)-block, and minus the identity on the (2, l)-block. Indeed, with 
respect to the symplectic splitting of Proposition |3.10| the basis 

ei = ( 1 1 ) , e2 = ( 1 1 ) , 

63 = ( 1 -1 ) , 64 = -L ( _i 1 ) , 

accomplishes this, as is easily verified. With respect to this basis, write the lineariza- 
tion at 7« = as 



L = JD2i7(0) = 



a b 
c d 



where a, 6, c, d are S0(2)-equivariant 2x2 matrices, which can therefore be identified 
with complex numbers. Since L is infinitesimally symplectic, d — —a* and 6, c are 
symmetric, i.e. d = —a and 5, c G M, and the eigenvalues of L are 

X = i Im(a) ± ■\/— Im(a)2 — (ac? — be) — i Im(a) ± \/Re(a)^ + be. 

So L is spectrally stable if the discriminant D = Re(a)^ + be is negative, unstable if 
D is positive, and an eigenvalue collision occurs if D = 0. Since D = when 1/ = 
and D < for 1/ < 0, ly ^ 0, D changes sign under the transversality condition 



dD 



7^0, (3.3) 

i/=0 
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whereupon the equilibrium becomes spectrally unstable for v > 0, v ~ 0, and under- 
goes a Hamiltonian Hopf bifurcation a,t v — see, e.g. [T7j- □ 

4 Application to the KirchhofF model 

In this section we use the theory from Section |3] to analyze the stability of vertically 
falling, spinning relative equilibria of a neutrally buoyant submerged axisymmetric 
body. In the Kirchhoff approximation (see e.g. [5J [7]), this is a Lagrangian system 
with configuration space SE(3), such that the position and orientation of the body 
in configuration {A,b) € SE(3) are obtained relative to a reference body, by the 
Euclidean transformation x '—^ Ax + b. The reference body is chosen with a vertical 
axis of symmetry, so that in the configuration (A, b) , the submerged body has axis of 
symmetry in the direction ^63. The Lagrangian is 

L{A, b, n, v) = ^n*in - mles ■nxv + ^v^Mv + mgl{e3 ■ Ae^). (4.1) 

Here (fi, v) S TSE(3) are the body-referenced angular and translational velocities, 
obtained by left translation, i.e. {^l,v) = {A,b)~^{A,b), or 

n = A-^A, V = A-^b. (4.2) 

The added inertia matrix I and added mass matrix M are the 3x3 diagonal matrices 

Ml 
M = \ Ml 
A/3 

and are determined from the shape and mass distribution of the body [6j . The reference 
body is such that the center of buoyancy is at the origin and the center of mass is at 
distance / below the center of buoyancy, so the vehicle is bottom heavy when I > 
and top heavy when I < 0. We will assume, as in the physical situation, that /, M , 
m and / satisfy 

Ji > 0, /3 > 0, Ml > 0, M3 > 0, /iMi - m^f > 0. (4.3) 




The Lagrangian is positive definite quadratic in under assumptions (4.3). 
The system admits the symmetries of 

(a) spatial isotropy: the left action of the subgroup SE(2) x R, corresponding to 
translations in any direction and rotations about the vertical; and 

(b) material symmetry: the action of the subgroup SO (2) by multiplication of the 
inverse on the right, corresponding to rotating the body about its axis. 

Thus the system has the symmetry Q considered in Section [3] The point x^ defined 

by 

Xe- A^id, 6 = 0, ri^^ea, v = ^e^, 

J 3 M3 
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is an axisymmetric relative equilibrium, corresponding to the motion where the vehicle 
spins at angular velocity Se/Ia and angular momentum Se about its (vertical) sym- 
metry axis, and translates along that axis with vertical velocity P^/M^ and vertical 
momentum P^. We will apply the general stability theory developed in Section [3] to 
this family of relative equilibria. 

4.1 Hamiltonian formulation of the Kirchhoff model 

We briefly review the derivation of the Hamiltonian and the equation of motion for 
the Kirchhoff model. See [6] for more details. The conjugate variables 11, P to A^b are 



— + mles XV, P = — 



= Mv — niles x fi, 

(A,6) = (id,0) 



and these equations can be inverted to give 



_ MiHi + mlP2 _ M1U2 - mlPi „ - 5^ 

_ hPi- min2 _ J1P2 + mffli P3_ 



The Hamiltonian is 



H{A, b, H, P) = {n, H) + {v, P)-L 



2{hMi - m^P) ^ ' 2h ^ 2(/iMi - m2p) 

^ -Pi + . ..""^^.,. (111^2 - n^Pi) - mglT, 



2M3 hMi - m^P 



where T = A e^. The equations of motion are 



(4.4) 



dt ^' dt ' dt 



This system is Poisson with bracket 



H^ 


pA 


pA 


pA 








pA 









(4.6) 



The momentum map J = (J**, J°, J^), c.f. (jOl, is given by 



J'^(A, b, H, P) = {AH + bx AP) ■ eg, 

J"(A,6,H,P) = ^P, (4.7) 
j\A,b,U,P) = -H-ea. 
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4.2 Reduction by x S0(2)r 

As shown in Section [3] the stability of the relative equilibrium Xe can be established 
by a study of the stability of the equilibrium of (2.7) at the parameter value v = 
(/Lig,/^g) = (PeB^jSe)- The system (2.7l is the Marsden-Weinstein reduction (see [8 ) 



of the full Hamiltonian system by the group M"^ x S0(2)fl;, at the momentum level 
(/^g,/^g), i.e. (2.7l is obtained from the full Hamiltonian system (2.11 by fixing the 



values of the conserved quantities J° and J to their respective values and fif, 
and eliminating variables along the group M.^ x S0(2)fl. For the Kirchhoff model, 
equations (4.5 1 cannot be directly used because they are the result of a reduction by 
the group S0(2)i x M^. In this section we compute the reduced spaces near Xe, in 



a way that relates to the variables 11, P, F of (4.5 1, by comoving with the action of 
S0(2)«. 

Given qi,q2, define q — qiBi + 9262 and 



A 



91:92 



id +x + fx' 



where / 



1 



x— — 63 X q 



(4.8) 



is orthogonal because / satisfies 



where x is defined in (2.2 1. The matrix in (4.1 

^/^ — 2/+ 1 = and / is the solution of this quadratic that is smooth 

it 



-63 X g is so that A^^^^^^e^ = qiBi + 9262 + 



if 



the equation — 2/ 

at q = 0. The choice x = 

The map {qi,q2,0) i— > Ag^^^^ exp(— 63^) coordinatizes S0(3) near the identity and 
provides an angle 0{A), which may be used to define the coordinates (a, 0, qi, 92) by 

{A,b) = {a,9) ■ {Aq^^q,_,mlM-\) = (A^^^,, exp(-6'e3), mZM-^g + a). 

The choice of b is such that the two sections a — 0,9 — and 91 = (72 = are 
orthogonal (at their intersection point) in the kinetic energy metric. This is an example 
of a general procedure advocated in [16]. The group x S0(2)fl acts by addition in 
a, 9, and 



F = exp(6'e3)(giei + 17262) + F3e3 



(4.9) 



so (91,(72) comoves with (Fi,F2). 

Let the conjugate momenta to a,9,qi,q2 be ly'^ , ,pi,p2. The coordinates a and 
9 are cyclic and so i^", i^^ are conserved. Thus the x S0(2)ij reduction results in a 
Hamiltonian in terms of canonical coordinates qi,Pi, which is parametrized by v°' . 
This is computed in the following proposition. 

Proposition 4.1. The Hamiltonian of the M"^ x S0(2)ij; reduction of the Kirchhoff 
model (|4.4| in the canonical coordinates qi,Pi is 



Hiqi,q2,Pi,P2) 



Mr 



2(/iMi - mH 

+ '-(^-- 

2 VM3 (/iMi -m2p) 

ml 



(pi+Pl - {12P1 - qiP2Y) 



(4.10) 



hMi - m^P 
mglT-i, 



{{qiP2- q2Pi){v'" ) + v^pi~vIp2) 
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wht 



Pi^ P2- fqiv 



ml 



P2 



-pi - fq2f 



ml 



(4.11) 



r = ^qi^^^^'s = 9iei +9262 -Tgea, ^3 = \Jl-ql-ql 



Proof. The x S0(2)ij reduced Hamiltonian can be computed by computing pi in 
terms of 11 and (71,(72 and eliminating P and li^ using (see ( |4.7[ )) P = A*;/" and 
Ha = — z/^. Since the Hamiltonian is invariant, it does not depend on a or 9, so we set 
a = and 9 = Q. By definition of canonical coordinates. 



Pi 



A 



tdA 
dqt 



n 



A' 



db 

dqt 



where the superscript ^ recovers a vector from an antisymmetric matrix by inverting 
the ^ operation defined by (2.2 1. Then 



At ^^gi,g2 

91,92 Qg^ 



= -62 



OA, 



91,92 



from which 



Pqiq2 



= 61- 



/-l 



63 X gr + Jq2 63 



Pi ^ 

P2=[l- 



/-I 



-1 



/-I 
/-I 



63 X (7 - /(7i 63 



n2 - fq2l^'^ 



/-I 



Hi 



Pqiq2 
/-I 



n2 + fqii^" 



ml 



These linear equations may be solved for 111,112, with the result 

^2 



Hi 

n2 



1 - fqi fqiq2 
fqiq2 1 - fql 



pi 

P2 



(4.12) 



where pi are defined by (4.11). The matrix in (4.12) is the same as the upper 2x2 
submatrix of —e^Ai „ 63, which is otherwise sparse, so setting p — piei +P262, 



and 



n = -63 X {Al^, Je, X p)) - z.«e3, (4.13) 



n? + n2 = ||-e3 X (yl*^_^^(e3 xp))|| = pj + pl - {q2Pi - qip2f , 

HiPz - n2Pi = 63 • (n X A\^^^.y) = [q^p^ - q2pi){v"- ■ r) + v^Pi - vlP2, 

^1-^-^^9263)^ 



The Hamiltonian (4.4) is obtained from (4.10) after substitution these and deleting 
inessential constants. □ 



17 



As already noted, the case of v"" vertical is important for the stability and KAM 
analysis, and we specialize to this now. 



Proposition 4.2. At vertical momentum 
system (4.4 1 is 

F, 



0, the reduced Hamiltonian 



H{qi,q2,pi,P2) = -^{{Pi + FiT^qi + v'^ fq2f - {qiPi + 92^2)^ 



+ {P2 + FiT3q2 - ly^fqi) 
+ ^hr + rngl (f- 



(4.14) 



-FpF^qW 



wht 



F„ 



Ml 



hMi -m?P' 



F, = mgl - {v^f 



1 

M3 



1 
Ml 



F 



Ml 



(4.15) 



Proof. Using (4.11), the Hamiltonian (4.101 can be written in the canonical coordi- 
nates (Ji, (?2jPijP2 as follows: 



H{qi,q2,Pi,P2) 



Ml 



2(/iMi - m^P, 

1 / J h 

2 \M3 (/iMi -m2p) 

ml 



Pi+pI~ ('72P1 - qihY) 



hMi ~ to2P 

Ml I 
2(/iMi-m2/2) I 

(z.f)2 / 1 



{qiP2 ~ q2Pi){i^3^3) - mglT^ 



Pi 



mlvii 
h 



12^:. 



mlvn 



P2 



Ml 

m2^2(i - ||g||2) 

M3 ^ (/iMi - 7712^2) ^ Mi(/iMi - m2/2) 



- (92P1 - 91^2)^ 

q\\'^ - mglTs 



Ml 

2(/iMi - m2/2) 



Pi + ^^91^3 + i^^92/ ) - {qiPi + q2P2f 



P2 



mlv^ 



-q2^z - V qif 



a\2 



+ - 2mglf - {v^) 



m^P\\q\\ 



M3 Ml Mi(/iMi - m2p) 



where a constant has been deleted at the final equality. Substitution of (4.151 then 
gives the Hamiltonian (4.14). □ 



Remark 4.3. In accord with Proposition 2.3 the Hamiltonian (4.14) admits an addi- 
tional SO (2) symmetry because v'^'' is vertical, which we can take as the diagonal ac- 
tion of SO (2) on ((71, 52), (pi,P2), with the standard conserved momentum qiP2 — q2Pi- 
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From (4.121, 



n • r = qipi + q2P2 

= qi{p2 ~ fqiv'^) + 92 (-pi - Iq2i''^) + n3r3 

= qiP2 - q2Pl - fy%ql + ql) - ^"(1 - \\qf f) 

= {qiP2 - q2Pi) - 1^^- 

Consequently, the additional conserved momentum (conserved for vertical v only) is 
equivalent to the subcasimir 11 • F of [7] . 

4.3 Energy-momentum confinement, spectral stability 
Theorem 4.4. The axisymmetric relative equilibrium Xe, with linear momentum 



and spin momentum Se, of the KirchhofJ model (4.51 is Q -stable by energy-momentum 
confinement if 



mgl > 



and spectrally stable if 

mgl > 



1 

Ms 



1 



1 

Ma 



1 



pf 



Ml 



4(/iMi - to2/2) 



Si 



(4.16) 



(4.17) 



Proof. By Proposition |3.1| it is sufficient to consider the equilibrium g = 0,p = 
of the M.^ X S0(2)i^ reduced KirchhofF Hamiltonian at the parameter values — Pe 
and — Se- The relevant Hessian may be obtained from (4.141 by deleting terms of 
higher than order two, resulting in 



H 



\p + Fq\' 



F, 



(4.18) 



where, for convenience, we use the notations 



qi+iq2, P^Pi + ip2, 



F^F, 



:Se 



By assumption Fp > 0, so the Hamiltionian (4.181 is positive definite when _F|j > 0, a 
condition which, in view of (4.151, is equivalent to (4.161. 



To show (|4.17|, note that the linearized equations at q — 0,p — are obtained 

's e 

dg2 



from (4.181 by Hamilton's equations 
dqi 



dq 
dt 
dp 
dt 



dt ^'~dt 



dpi 
dt 



dH 
dpi 



.dH 

dp2 



dp2 



dH 
dqi 



.dH 

' dq2 



FpFq + Fpp, 

-iFp\F\^ + Fg)q^FpFp. 
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Due to S0(2)-equivariance, computing the eigenvalues of the linearization reduces to 
solving the equation 



= det 
This has roots 



F„F-\ 



-{Fp\F\^ + Fq) -FpF^X 



= X^ + iFpSeX + F^Fa. 



\=\ [~iFpS, ± ^-F^Sl - 4FpF,) , 



so spectral stability holds if and only if 
+ 4FpFq = AFp (mgl 



1 1 

M3 " Ml" 



Ml 



4(/iMi - m^P) 



(4.19) 



> 0, (4.20) 



which is equivalent to (4.171. 



□ 



Remark 4.5. Energy-momentum confinement under momentum-preserving perturba- 
tions, and spectral stabihty, fail coincidentally as predicted in general by Proposi- 
tion 3.11 and as observed in this example by [7]. We can check this by determining 
the conditions such that H + Fp\{qip2 — q2Pi) is positive definite for some A € M, 
where H is the Hessian in (4.18). One computes 



H + X{qiP2 - q2Pi) = 



^{Fp\p + Fq\^ + 2X{q,p2 - <72Pl)) + 

^jFpip + Fq) + zXq\' + ^ {FpF, + XFpS, - X') \q\' , 



(4.21) 



where we have used 

\Fpip + Fq) + zAgp = F^\p + Fq\^ + X^\q\^ + 21m{Fp{p + Fq){Xq)-) 

= F^\p+Fqf + X'lql^ + 2XFp{{qiP2 - q2Pl) + |-z|' Imi^) 



In (4.21 ), the coefficient of \q\ is positive for some A if and only if FpFq + XFpSe — X 
has a real root, i.e. if and only if 

FlSl - 4{FpFq){-\) = FlSl + AFpF, > 0, 



which is exactly the condition for spectral stability as in (4.20). 



4.4 KAM stability 

In this section we verify that (at vertical j^°) the equilibrium q = 0,p = of the 
X S0(2) reduced Hamiltonian ( |4.10 1 satisfies the twist condition of the Arnold 
stability theorem (see [9]) in the gap between (4.16) and (4.17): 
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Theorem 4.6. The axisymmetric relative equilibrium Xe of the Kirchhoff model (4.5 1 
(with linear momentum and spin momentum Se) is stable if 



pz 



1 



1 

Ml 



> mgl > 



1 

M~3 



1 

Ml 



Ml 



4(/iMi - m^P 



For the proof we use Arnold's Theorem, see e.g. Chapter IX. E]. This requires us 
to bring the Hamiltonian into the normal form 



H = Ho 



2N 



(4.22) 



where 



(a) H2k, 1 < £ 27V are homogeneous degree k in the polynomials Ti = ^{qf +pl) 
and J2 = 

(b) is at least order 2N + 1 in qi,q2,pi,p2', 

(c) H2 = uJiXi — LU2I2, where lui,lu2 =/= 0. 

Theorem (Arnold Stability Theorem). The origin is stable for the Hamiltonian sys- 



tem (4.221, provided for some k, 2 < k < N , the twist condition I?2fe = -ff2fc(w2,wi) 7^ 



is satisfied. 



Proof of Theorem 4-6 Necessarily u ~ Se ^ 0, or else there is no gap. By choice of 
time scale, we can assume that Se = 1. Also, Fq < in view of (4.151 and (4.171. The 
proof requires four steps. 

1. Consolidate parameters. The conditions D2k 7^ of the Arnold Stability Theorem 
are algebraically complicated in the parameters. The number of these conditions that 
is expected to be required is one more than the number of free parameters in the 
Hamiltonian, so as to have the reasonable expectation that there will be no value 
of the free parameters for which all conditions vanish. To facilitate the symbolic 
computation of the required normal form, so that the computation to be tractable, it 
is necessary to reduce as much as possible the number of parameters, and to cast the 
remaining free parameters so that they appear as simple powers. 



First, the parameter Fi of the Hamiltonian (4.141 may be eliminated by the sub- 
stitution 



Pi = ui 



Fm 



P2 = U2 



Fiq2 



This is a symplectic transformation because it is a shift by the exact form 



FidV^ 



Fiqi 



dqi 



Fiq2 



dq2 



By a simple direct computation. 



{pi + FiT^qi + fq2Y - {qiPi + 92^2)^ + (P2 + FiT:iq2 - fqiY 
= (^1 + fq2f ~ {qmi + q2U2f + {U2 - fqi? - 
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so (4.101 is transformed to 



where, evidently, Fi is absent. 



The frequencies of the hnearization of g = 0,p = are related to Fp, Fq by (4.191: 



fi = liFp + ^F^ + AFpFq), f,= l(Fp^^F^+AFpFq), 



and these can be inverted to give 

Fp — fi + f2, Fq 



/1/2 

fl+f2' 



(4.23) 



Note that /i, /2 > and /i - = + AFpFq > so that /i > /2 > 0. Replacing 

mgZ with — /i, and deleting an inessential constant, it is sufficient to consider the 
Hamiltonian 



fi + /2 

ff(gi,92,"i,M2) = — 7, (("1 + fl2)^ - ili'^i + 92^2)^ + {u2 - fqiT) 



(4.24) 



under the assumptions that 



/l > /2 > 0. 



The number of free parameters is 2, because an inessential multiplier of the Hamilto- 
nian can reduce the triple /i, /2, /i by one. Thus, three twist conditions are expected 
to be required. 

2. Linear normal form. The linearization at the origin with respect to the coordinates 

(gi,g2,ui,U2) is 



'4(/l+/2) 





\{h + h) h + h 

/1 + /2 

i(/l+/2) 



/ 2U1 + J1J ./H-/2 

\{h + h) 

(/i;:/2)! Q 

(/1-/2) ' 

4(/: 



and four linearly independent eigenvectors are 



ei - ( 2(/i + /2) /2 - /i ) 

62 = ( 2(/i + h) /i - /2 ) 

63 = ( 2(/i + /2) /2 - /l ) 
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64 - ( 2(/i +h) /l - /2 ) 

with respect to which the symplectic form is canonical with multiplier 4(/i — /2)(/i + 
/2), which can be ignored. The coordinates {Qi, Pi,Q2, P2) defined by 



( qi (?2 Ui U2 



normalizes the quadratic part of H to 

2 



)iei + Pi 62 + (5263 + -P2e4 



where 



= 4(/i - /2)(/l + /2)/l, ^2 = 4(/i - /2)(/l + /2)/2. 



The S0(2) action for vertical v explained in Remark 4.3 is anticlockwise rotations on 
{Qi,Pi) and clockwise rotations on ((32,^2)- The momentum of this action can be 
taken as 



-X1+X2. 



(4.25) 



3. Nonlinear normal form. Since the Hamiltonian H is invariant under an S0(2) 
symmetry, we choose to seek normalizations within the corresponding class of invariant 
functions. Convenient invariants are 



1 



W3 



-{Q1Q2-P1P2), 



W2 
W4 



-{Q1P2 + Q2P1), 



and these satisfy the relation 

2WlU>2 — Ulg — W4 = 0. 

The Poisson bracket of the Wi is closed and the matrix {wi, Wj} is 



(4.26) 



{Wi,Wj} 








W4 


-W3 










-W3 









-Wi - W2 


\ W3 


W3 


Wi + W2 






(4.27) 



Since H is invariant it can be expanded in a Taylor series in wi, ti;2 7 W3, W4 as 
H = H2 + Hi + --- + H2N+H'^ 

where 

H2 = LUlIl - UJ2I2 = UJlWl - UJ2W2 
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and the are degree i homogeneous polynomials in wi, W2, ws, W4. The normal- 
ization process is formally the same in the space of invariants as it is in the space 
of functions of qi,pi,q2,P2- See [9j Chapter VII. C]. The normal form is achieved at 
degree k > 2 successively in k. After normalizing the terms of degree smaller than 
or equal to k, at the beginning of step k, we start with the Hamiltonian H2k from 
the previous step. The Hamiltonian i?2fc is in normal form up to order k, i.e. it is 
the sum of homogeneous degree i < k polynomials in wi , W2 and homogeneous degree 
i > k polynomials in wi,W2,W3,Wi. We now look for a homogeneous polynomial G 
of degree k + 1 in wijW2,W3,W4 such that the time 1 flow of G transforms H into a 
function 

H = H + {G,H}+^{G,{G,H}} + --- , 

which is in normal form at degree k+l. The degree k+1 term of is i?2(fc+i)+{G, H2} 
because 

deg{G, H2^} - degDG + degDi?2, + l^k + {i-l) + l^k + i, 
so deg{G, H2i} = fc + 1 only for i = 1. The coefficients of G are adjusted so that 

H2k=H2k + {G,H2} (4.28) 
is a function only of wi and W2, after all powers of W4 greater than 1 are eliminated 



using (4.26 1. The Poisson brackets here are computed using (4.27). Then H is replaced 
by H and the computation proceeds to step fc + 1. 

At step k, there always exists a G such that H2k is normalized. Indeed, {-,^^2} 
maps the finite dimensional vector space Pfc+i of homogeneous degree fc+l polynomials 
to itself. The space Pfe+i can be regarded as an inner product space by 

fg 

B12 X -B34 

where i?i2,-B34 are the unit disks in the wi,W2 and W3,W4 variables. The linear map 
{ • , H2 } is antisymmetric with respect to this inner product because 



{f,{9,H2}) + {{f,H2},9)= [ {fg,H2} 

B12 X -B34 



. .f difg) ^ d{fg)\ 
[uji - u;2j -Wi— h W3— 

B12XS34 V ^■'^3 OWi J 



= -(cji-cj2) / d> fgiwsdws + Widwi 



= 

because the vector field in the line integral is always orthogonal to its path. Conse- 



quently, the image of {•,772} is orthogonal to ker{-, 7J2} and G in (4.28 1 can be chosen 



so that B.2k € ker{-, i72}. Moreover, / G Pfc+i is in ker{-, H2} if and only if 
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from which / = f{wi,'W2, + w\), since loi ^ lo2- But in view of the relation (4.261, 
this implies that / is a function of wi, W2, as required. We note that only the nonreso- 
nance condition wi ^ UJ2 is required and the normalization could in theory be achieved 
at any order. 

The normal form cannot be hand-computed, but can be computed using a symbolic 
manipulator. The result up to and including iJg (the notation 1 2 means the 
preceeding fragment with wi and W2 exchanged) is 



Hi = 8(/i + f2f[A2owl + ]^A^iWiW2 + 1 ^ 



32(/i + 12? 

(/l - h? 

64(/l+/2)V4 .,,4 , . .3 , 1 , 2 2 



H. - 



Hs = - (^J _ J Y ("^40^1 + AiiW-yW2 + -^42^1^2 + 1 2 



A20=(-/l-/2)Ai+(/l + /2)/2, 

A21 = (-4/1-4/2)^ + 8/1/2, 

A30 = 2(/l + /2)V' - 2(/l + /2)(3/l + /2)m/2 + 4(/l + /2)/l/2^ 

A31 = 15(/i + /2) V' - (/i + /2)(5/i' + 44/1/2 + 11/2')/^ 
+ (5/2 + 38/1/2 + 1 7/1) /1/2, 

^32 = 15(/l + /2) V' - (/l + /2)(ll/f + 44/1/2 + 5/|)/i 

+ (17/2 + 38/1/2 + 5/|)/i/2, 
A40 = -16(/i + /2)V' + 2(/i2 + 36/1/2 + ll/|)(/i + /2)V' - 2(/i + /2)(2/f 
+ 55/2/2 + 36/1/2 + 3/3)^/2 

+ 2(/l + /2)(/i2 + 26/1/2 + 5/22)/l/22, 

A41 = -182(/i + /2)V' + 26(3/2 + 31/1/2 + %fl){h + /2) V 

- 2(13/i + /2)(/i + /2)(8/i2 + 49/1/2 + 21/22)m/2 

+ 2(65/3 + 367/2/2 + 267/i/2 + 29fi)hfl 
A42 = -354(/i + /2) V' + 3(95/2 ^ ^i<ifj^ + 95/22)(/i + f2fpL^ ~ 3(/i + /2) 

(9/4 + 274/3/2 + 850/2/22 + 274/1/I + %fl)ii + 3(9/4 ^ 204/3/2 

+ 518/2/1 + 204/1/3 + 9/4)/i/2. 

^. Twist condition The twist condition in the Arnold stability theorem is D2k — 
H2k{^2i^i) 7^ for some k >2. From the normal form, the first three Z?2fc are 

Di ^ 128(/i - /2)2(/i + /2)'(-(/l + /2)(/l' + 4/1/2 + /2')M 

+ (/? + 10/1/2 + /|)/l/2), 
D, = -2048(/i + /2)^((2/i2 + 13/1/2 + 2/22)(/l + ./2)2m2 - (/i + /2)(ll/i2 

+ 46/1/2 + 11/2')m/i/2 + (S/r + 50/1/2 + 9.fl)flf2)ifl - /2), 
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= 16384(/i 



/2)"(-2(8A4 + 91/3/2 



91 /i /I 



f2')(/l+/2)V 



177/f/l 

(2/f + 150/f /2 + III3/1V2 + 1970/3/1 + 1113/2/1 
I5O/1/I + 2/|)(/i + /2)V' - (/i + /2)(4/f + 345/f /2 + 2226/iVI 
3850/f/l + 2226/2/1 + 345/i/| + 4.f^)tifif2 + (2/f + 211/^/2 
1466/f /I + 2642/3/1 + I466AV2 + 211/i/| + 2/|)/2/|). 



We assume Dq, Dg vanish and argue by contradiction. The expressions are homo- 
geneous in /i, /2, so one can assume /2 = 1- Since D4 is hnear in if the coefficient of 
H vanishes, then the assumption = imphes the constant term in fi also vanishes, 
i.e. 

-(/i + l)(l + 4/i + /2) = 0, and /i(/2 + 10/i + l) = 0. 

But these two polynomials in /i cannot simultaneously vanish because their gcd is 1, 
as can be computed with the Euclidean algorithm. So the coefficient of fi in D4 does 
not vanish. Solving D4 for /i and substituting into Dq and Dg results in 



= 10240(/i + iffHh + 5)(5/i + l)(/i - 1)3 
= -245760(/i + l)"/i4(21/4 + 172/3 + 334/f 



172/1 + 21) (/i-l) 



But these cannot simultaneously vanish because their gcd is 10240(/i + l)^/f (/i — 1)3, 
which is nonzero in view of /i > /2 = 1. □ 



Remark 4.7. The computer generated normal form used in Theorem |4.6| requires an 
independent check, and one possible check follows. The reduction by SO (2) of the 



S0(2)-invariant Hamiltonian (4.241 gives a one parameter family of 2-dimensional 



Hamiltonian systems parametrized by the S0(2) angular momentum. All trajectories 
of these reduced Hamiltonian systems near q = 0,u = are periodic. The periods 
of these orbits depend only on the energy and the S0(2) momentum, and may be 
numerically computed. On the other hand, the successive normal form Hamiltonians 



nf . 



H2 + H4 + Hg 



(4.29) 



iJg — H2 + H4 + Hq + Hs 



are functions of wi and 'W2- Using the Poisson bracket (4.27), the differential equations 
for Wi are 



dwi 



dW2 



= 0, 



dW3 

dt 



dwi 



dim 

dw2 



W4, 



dw4 
lit' 



dwi dw2 



W3. 



These differential equations are the S0(2) reduction of the Hamiltonian system (4.241 
and the orbits have periods 



■nf 



dwi dw2 



(4.30) 
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mm 


2i = 4 




2i = 6 




2i = 8 




1 


.9995410553688 


.99954076 


4.0 


.99954105570 


6.0 


.9995410553684 


8.2 


2 


.9981715477722 


.99816682 


4.0 


.99817156906 


6.0 


.9981715476563 


8.0 


4 


.9928005386538 


.99272712 


3.3 


.99280185244 


5.8 


.9928005101406 


7.8 


8 


.9728669887964 


.97182458 


3.5 


.97294031478 


5.4 


.9728607222126 


7.2 


16 


.9108473444693 


.89972618 




.91385268217 




.9098792303872 





Table 1: Comparison, at the same energy-momentum values, of numerically computed 
periods, and periods computed from the norm,al forms i?2f . The correct normal form 
is indicated because the values of r in the columns are nearly the corresponding values 
of2i. 



which also depend only on the energy and the SO (2) momentum. For energy and 
momentum obtained from initial conditions eq^ , ep^, the periods of -ff2i^, computed di- 



rectly from (4.30 1, can be compared with numerically computed periods of the Hamil- 
tonian (4.241. If the normal form is correct, then the difference of these must fall as 
e^*, because the period from ifj/ is accurate to order i — 1 in w, and therefore its 
difference falls as order i, which is order 2i in e. If the normal form is not correct, 
then it is not likely that the anharmonic periods it predicts will agree at the proper 
order with the periods of its non-normalized precursor. 
We have implemented this check for the data 



Vb 1 

-y, /2 = 1, A^=2, 1 
6 = 0.0001, 0.0002, 0.0004, 0.0008, 0.0016 



/: 



(0.5,1.0), / = (-0.75,0), 



The choice of /i is convenient because then 4(/i — /2)(/i + /2) = 1 and the symplectic 
multiplier of Theorem |4.6| does not have to be accounted for. In Table [l] are the 
numerically computed periods T{e), and the periods T'^/(e) computed using the normal 
forms, as ratios of the period 



Tn 



2n 



2tt 



/1-/2 



2 



1 



of the linearization. Also tabulated are the values called r2i^e, which are the base 
two logarithms of the successive (as e is successively doubled) ratios of the differences 
between the periods r^/(e) and T(e), i.e. 



r2i 



log2 



T2f(2e)-T(2e) 



T-^{e)-T{e) 



In the column corresponding to H2i, there is agreement between r2i,e, which is the 
observed order that the difference falls as, and 2i, which is the predicted order. We 
conclude that there is numerical evidence that the normal forms iJj/ are correct, 
because they predict the correct anharmonic period with an error of order 2i. 
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5 Numerics 



In this section we simulate small dissipation directly in the M'^ x S0(2)ij reduction of 



the KirchhofF model, as computed in Proposition 4.1 We verify that, after addition 
of small dissipation, 

(a) stability is maintained for the axisymmetric relative equilibria Xe which are in 



the EM-region determined by (4.161; 



(b) stability is destroyed for the axisymmetric relative equilibria Xe which are in the 
gap determined by (4.161 and (4.17). 



This confirms that the stability of the system in the gap is qualitatively less robust 
than its stability in the EM-region, and therefore confirms that energy-momentum 
confinement fails in the gap, with physically significant consequences. The approach of 
simulating the xS0(2)fl reduced system, rather than the full KirchhofF model (4.1 1, 
is consistent with our development in Sections |2] and [3] and is convenient for the 
simulation of energy dissipation which preserves the M.^ x 80(2)^ momentum and 



symmetry (see Remark 3.7) 



5.1 Splitting method 

The Hamiltonian ( |4.10 1 may be split as 

H = Hq + Hi + H2 + 



where 



Ho 
Hi 

H2 

H3 



Hiqi,q2,0,0), 



pI - 2q2P2qiPi - qfpl - qjpl), 



FpFi{qipi + 92P2)Y 1 
Fpiy^{p2qi -~Piq2) 



1 + 



9? 



91 



Only Hq depends on Va, and all but Hq are SO (2) invariant under the diagonal action 
of S0(2) on (gi, (72), (pi,P2), whereas Hq has this invariance only for vertical v (see 
Remark [4!3|. 

There are exact formulas for the flows of each Hi, as follows: 

(a) Hq is a function of q only and hence its flow is p p — tV Hq. 

(b) The flow of Hi corresponds to the motion of a particle of mass l/Fp freely 
moving on a sphere of radius 1. Indeed, for the top half of the sphere with 
coordinates x = qiei + 5261 + ^1 — 9^ — g| 63, the kinetic energy metric of the 
particle motion, and its inverse, are respectively 



1 



FpT^ 



9192 

1-9? 



1 



9? 

-9192 1 



-9192 

4 
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the second of which is twice the matrix corresponding to the quadratic form Hi . 

(c) Since H2 is Hnear in p, and pi,P2 are separated, the differential equations for 
qi,q2 close. As is easily verified, qi/q2 is conserved for these equations and 
so qi,q2 evolve along radial lines. The radial equation can be separated and 
integrated, and the evolution of pi , p2 follows that from conservation of energy 
and momentum. 

(d) As is easily verified, qf + is conserved for the flow of H^. Using this, the 
differential equations corresponding to become linear and may be integrated. 

For a dissipation vector field, one could use a positive multiple of the negative 
gradient in pi,p2 of the reduced Hamiltonian (4.101, i.e. the vector field 



(1 - + q,q2P2 - -^q, " j ^ 



+ \qxqiVx - (1 - 92)P2 - ^^92 + ^^<Zl j 



dV2 



(5.1) 



This does not preserve the momentum q\P2 — P\q2 which occurs for vertical v (see 
Remark 4.3 1. But the direction of (5.1) relative to the kerDi? hyperplane (for <7i,g2 
constant) is unaltered by projection: any projection of (5.1) onto any subspace will 
still dissipate energy. So we project onto the line tangent to the level set of the function 
qiVi — Viqi, i-e. onto the vector (51,92)- Since the inner product of (5.1) and (91,92) 
is 

(-1 + 9? + 9^)9ipi + (-1 + 9? + (i)q2V2 - "^^Jy^(9? + g|) 

Ml 

= -^{qiVi + <Z2P2) - + ll)^ 

we use for the dissipative perturbation vector field ei?, where e is small and positive, 
and 

d 

R = {-{qiPi + q2P2)T3 - F;(9i + 92)) qi^ 

Q (5.2) 

+ {-{qiPi + 92P2)r3 - Fi{ql + ql)) 92^- 

The differential equations corresponding to this vector field are linear and easily solved 
exactly. 

The full system, dissipation included, is given by the vector field 



l.S 



dH, d dH, d 



dpi dq^ dqi dpi 



(5.3) 



Letting the flows of the Hamiltonians 77^, < s < 3, be and the flow of the 

dissipation ei? be Gt, the concatenation 

P^O I7'^2 I7'^2 

'J"At/2-f^ At/2 At/2 At/2-^ At At/2-^ At/2-^ At/2'-^At/2 
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is a second order one step method, which we use to simulate the system. At vertical 
this method preserves the conserved quantity qip2 — q2Pi because every one of its 
steps does that separately. The method is symplectic if e = 0. 



5.2 Simulations 

The simulations that we report in this section are done at parameter values 

2 (5.4) 
i.e. units are such that g — I. The initial conditions are chosen so that 

92 = 0, pi= 0, P2 = 0, = 0. 

The perturbation from the equilibrium is achieved by choosing small nonzero qi and 



From (4.4 1, energy-momentum confinement occurs for < Ci and spectral 
stability occurs for < C2, where 

M1M3 M1M3 / Ml , 

^^=Mr^™^^' ^ + 4(/iMi - m^P) 



For the parameter values (5.4 1, Ci = 1 and C2 = 4, so 



2 < Pe : spectral instability; 

1 < Pe < 2 : gap; spectral stability; 
< Pe < 1 : EM-region; stability. 

5.2.1 Stability 

On the left. Figure [l] shows a regular time sampling of the radius r = \/q\ + q2 vs. 
time. The top and bottom graphs correspond respectively to Pe = 1.5 (which is inside 
the gap) and Pe = .5 (which is inside the EM-region). The graphs have been obtained 
from 

e = .05, (7i==0.05/r, ul = Qm/k, fc = 1, 1.5, 2, 3, 4, 

e = .1, q^ = 0.05/r, = 0.01/fc, k = ^2, {V^f , (^/2)^ (^/2)^ (72)1°. 

and for 3 x 10^ and 1.2 x 10^ periods, respectively. Here, "period" means the smallest 
period of the normal modes of the linearization at the equilibrium q = p ^ Q, i.e. 271 /lo 



where uj is the largest of the imaginary parts of the eigenvalues computed from (4.191. 
The top and bottom of Figure [T] corresponds to the first and second lines of initial 
data above, respectively. The time step for the simulations was At « .04453. The 
simulations were stopped if r exceeded .5, corresponding to a vehicle configuration 
which is skewed 30 degrees to the vertical. 
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1.08 X 10*^ 



1.08 X 10*^ 




AH 



.004 



3.5 X 10^* 




3.5 X 10" 



Figure 1: Left top: instability in the gap. Left bottom: stability in the EM-region. 
Dissipation is added in both the top and bottom, and the energy decay is shown in the 
graphs at the right. 



At the top left, instability is indicated in the gap because large deviations from 
the equilibrium are observed for sufficiently long time. At the bottom left, stability is 
indicated in the EM-region because of the decreasing deviation from the equilibrium 
for decreasing initial conditions, irrespective of the elapsed simulation time. At the 
right, energy dissipation is demonstrated by plotting the energy minus the initial 
energy (Ai7) against time. 

5.2.2 Transitions 

On the left bottom. Figure |2] shows the maximum value of r over 1.2 x 10^ periods, 
for 



e = 0, gi= 0.0125, 



0.0025. 



At was adjusted to 40 time steps per period. On the right is the same except for 
e = .05. For the no-dissipation runs on the left, stability in the gap between Ci and 
C2 is indicated by the small maximum values for r. On the right, large deviations 
from the equilibrium occur over the entire gap between Ci and C2. 

The peak at Ci in the left graph indicates loss of stability at that transition, even 
in the purely Hamiltonian context. This is due to the presence of a zero eigenvalue 
in the linearization (see (4.191 when Fq = 0). At Ci, KAM stability as discussed in 
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Oii=it^- , , 1 0^>=^E^- , , 

12^^ 1 2 

Figure 2: Destruction of KAM stability of an axisymmetric relative equilibria in the 
gap between Ci and C2 as a result of adding small dissipation. Left: no dissipation, 
right: added dissipation. Top: plots the real parts of the spectrum of nearby persisting 
relative equilibria. The scale on the lower parts of the graphs differs from the scale of 
the upper parts so that the sign change at Ci is visible. The computation was done with 
the program MAPLE at 18 digits of accuracy. Bottom: large r indicates instability. 



Sections |3.1| and |4.4| is not present because it requires perturbation from an elliptic 
equilibrium but the equilibrium is not elliptic, it has a eigenvalue. 



The equilibria of the w-equation (5.3 1 aX q = p = Q and vertical momentum v 



persist to nearby nonvertical momenta. As further evidence of dissipation induced 
instability, we have computed the corresponding spectrum of the linearization of these 
equilibria, which were found numerically by Newton's method with start at (7 = p = 0. 
As shown in the top of Figure |2] at zero dissipation, we observe zero real parts in the 
spectrum throughout both the EM-region and the gap. For nonzero dissipation, the 
spectrum splits and has negative real parts in the EM-region and positive real parts 
in the gap. The real parts in the gap are small as compared to the real parts after the 
Hopf bifurcation, necessitating the two vertical scales in the Figure. 
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